Redshift drift in a pressure gradient cosmology 
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We derive the redshift drift formula for the inhomogeneous pressure spherically symmetric 
Stephani universes which are complementary to inhomogeneous density Lemaitre-Tolman-Bondi 
(LTB) models. We show that there is a clear difference between the redshift drift predictions for 
these two models. The Stephani models have positive drift values at small redshift and behave qual- 
itatively as the ACDM models while the drift for LTB models is always negative. This prediction 
can be tested in future telescopes such as E-ELT, TMT, GMT, and especially, in gravitational wave 
interferometers DECIGO/BBO. 
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In the context of dark energy problem there has been 
more interest in the non-Friedmannian models of the uni- 
verse which could explain the acceleration only due to 
inhomogeneity [ll, One of the strongest claims was 
that we are living in a spherically symmetric void of den- 
sity described by the Lemaitre-Tolman-Bondi (LTB) dust 
spheres model Q. However, there is a variety of inhomo- 
geneous models (for a review see e.g. Q) which have the 
advantage that they are exact solutions of the Einstein 
field equations and not the perturbations of the isotropic 
and homogeneous Friedmann models. One of the reasons 
to investigate the simplest models is their mathematical 
feasibility supported by the paradigm of the Copernican 
Principle which says that we do not live at the center 
of the Universe. However, observations are practically 
made from just one point in the Universe ("here and 
now" ) and extend only onto the unique past light cone of 
the observer on the Earth. It is clear that even the cos- 
mic microwave background (CMB) is observed from such 
a point. Apparently, its observations prove isotropy of 
the universe (isotropy with respect to observation point - 
the center of symmetry), but not necessarily homogene- 
ity (isotropy with respect to any point in the Universe). 
Then, the question is whether we should first start with 
model-independent observations of the past light cone 
and then make conclusions related to modeling of the uni- 
verse (cf. observational cosmology programme of Ref. [sj ) . 
In other words - homogeneity needs a check. Suppose 
that we have an inhomogeneous model of the Universe 
with the same number of parameters as a homogeneous 
dark energy model and they both fit observations very 
well. How could we differentiate between these two mod- 
els? 

The simplest inhomogeneous cosmological models are 
the spherically symmetric ones and these are the comple- 
mentary to each other: the inhomogeneous density g(t, r) 
(dust shells) LTB models and the inhomogeneous pres- 
sure p{t, r) (gradient of pressure shells) Stephani models. 



Apparently, for some reasons most of the cosmologists 
investigate the former, and only a few investigate the 
latter. In this paper we would like to investigate this 
complement of an LTB models - the spherically sym- 
metric Stephani model of pressure gradient which also 
possesses a totally spacetime inhomogeneous generaliza- 
tion d, 0] ■ Actually, these were the first inhomogeneous 
models ever compared with observational data from su- 
pernovae [1] and showed that they can be fitted. Despite 
the LTB models were theoretically explored much earlier, 
only later they were tested observationally @. 

For the sake of comparison we remind that the 
Lemaitre-Tolman-Bondi universe is the only spherically 
symmetric solution of Einstein equations for pressureless 
matter energy-momentum tensor T"*" = gu°'u^ [g is the 
energy density, is the 4-velocity vector) and no cos- 
mological constant which has a spatially dependent " cur- 
vature index" k{r). It has to obey some "regularity con- 
ditions" like: the existence of a regular center of symme- 
try, the orthogonality of hypersurfaces of constant time 
(of topology which implies the existence of a second 
center of symmetry) to a 4-velocity vector, and the avoid- 
ance of a 'shell-crossing" singularity (see e.g. [3]). Kine- 
matically, they are characterized by the non-vanishing 
expansion scalar O and shear tensor Gab- One of the pe- 
culiarities is that in LTB models the big-bang singularity 
is not necessarily instantaneous - different points start 
evolution in different moments of time. 

On the other hand, the inhomogeneous pressure 
Stephani model is the only spherically symmetric so- 
lution of Einstein equations for a perfect-fiuid energy- 
momentum tensor T"** — {g + p)u°-u^ + pg"^^ [p is the 
pressure, g""^ is the metric tensor) which is conformally 
flat (Weyl tensor vanishes) and embeddable in a 5- 
dimensional fiat space [1, [3|- A general model has no 
spacetime symmetries at all, but its 3-dimensional hy- 
perspaces of constant time are maximally symmetric like 
in the Friedmann universe. In order to be consistent with 
an LTB, here we consider only its spherically symmetric 
subcase which reads as (one uses a Friedmann-like time 
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coordinate ^) 



+ ^[dr^ + r^{de^ + sin^eV)] , (1) 



where 



V{t,r) = l + -k{ty 



(2) 



and (...)■ = d/dt. The function a{t) plays the role 
of a generalized scale factor, k{t) has the meaning of a 
time-dependent " curvature index" , and r is the radial co- 
ordinate. Kinematically, these models are characterized 
by the non-vanishing expansion scalar Q and the accel- 
eration vector 

The energy density and pressure are given by 



p{t,r) 
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V(t,r} 
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Vit,r) 
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(3) 



g{t)=We{t,r)g{tl4,) 



where we have assumed the velocity of light c = 1, G is 
the gravitational constant, and We{t, r) is an effective spa- 
tially dependent barotropic index. It is important that 
the metric ([T|) is written in terms of the isotropic coordi- 
nate which is related to a standard Friedmann coordinate 
f via the transformation f — 1 + {l/2)kr'^. The global 
topology of Stephani models is x R and they look like 
the de Sitter hyperboloid with specific deformations near 
the "neck circle" Q- The local topology of the constant 
time hypersurfaces (index k{t)) may change in time. In a 
standard de Sitter case one cuts the hyperboloid by either 
k^l (S^ topology), k = {R^ topology) or fc = -1 (iJ^ 
topology). Here we have a "3-in-l" case and the Uni- 
verse may either "open up" to become negatively curved 
or "close down" to become positively curved. In a gen- 
eral Stephani model which has no spacetime symmetry 
at all, the point reflecting an instantaneous ("only one 
hypersurface" ) center of symmetry moves around a de- 
formed hyperboloid. Similarly as in an LTB model, there 
are two antipodal centers of symmetry. In the Stephani 
models there exist instantaneous standard big-bang sin- 
gularities (a— >-0,p— ?>oo,p^ooas well as finite 
density (FD) singularities of pressure which appear at 
some particular value of the radial coordinate r 0, [Tlj . 
FD singularities at first glance resemble sudden future 
singularities (SFS) [l^ which appear in Friedmann mod- 
els with no equation of state to link the energy density 
and pressure. However, FD singularities occur as singu- 
larities in spatial coordinates rather than in time, which 
means that even at the present moment of the evolu- 
tion they may exist somewhere in the Universe. It has 
been shown in Ref.[20| that SFS may appear in inho- 
mogeneous Stephani universes, independently of the FD 



singularities. In Stephani models there is also the space- 
like H-boundary which divides each negative curva- 
ture k{t) < hypersurface onto the two sheets (the "far 
sheet" and the "near-sheet" @). It appears whenever 
the function V{t, r) in ([2]) is zero. On a H-boundary the 
Universe behaves asymptotically like de Sitter. As one 
can conclude from ([3|) and Q , there is no global equation 
of state, rather it changes from shell to shell. 

In Refs.0, two exact spherically symmetric 

Stephani models were found: the model I which fulfills 
the condition iV/a)" = 0, and the model II which fulfills 
the condition (k/a)' =0 (this reduces the factor in front 
of dt"^ in the metric ([T]) just to —1/V^. A subclass of 
the model II with k{t) = fia{t) {(3 = const., with the unit 
[/3] — Mpc~^) was found in Ref. To'] and it assumed that 
at the center of symmetry the standard barotropic equa- 
tion of state p{t) — wg{t) was fulfilled. This assumption 
gives that 
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git) 



{A — const.) 



(5) 



3 a3("'+i)(t) 
and allows to write a generalized Friedmann equation as 



d{t) 



A' 



a{t)J a3("'+i)(t) a{t) 
with the equation of state 



(6) 



Pit) = 





w + ^{w + l)a{t)r' 



g{t) = Wegit) 



(7) 



Similarly as in the Friedmann models, one can define 
the critical density as gcr{t) — {3/8TrG)[d{t)/a{t)]'^ , and 
the density parameter ri(t) = g{t)/ gcr{t)- After taking 
t — to, we have from ([5|) that 



1 = 



A^ 
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— Qo + ^inh , (8) 



and so 



f3 = aoH^ [no - 1) 



(9) 



The non-vanishing components of the 4-velocity and 
the 4-acceleration vectors are 



V ' 

The acceleration scalar is 



a 



V 



2^ ' 



(10) 



(11) 



and it does not depend on the time coordinate at all. 
Bearing in mind that the constant /3 is negative in our 
model (cf. formula ([9]) for ilo < 1), we have that the high- 
est pressure is at r = (the center of symmetry) , while 
the lower (negative) pressure regions are outside the cen- 
ter, so that the particles are accelerated away from the 
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center which is similar effect as the effect of the posi- 
tive cosmological constant in ACDM model. However, 
in ACDM the pressure is constant everywhere while in 
the spherically symmetric Stephani model it depends on 
the radial coordinate r. The components of the vector 
tangent to zero geodesic are [HI 



J.2 ' 







(12) 

where h — const., and the plus sign applies to a ray mov- 
ing away from the center of symmetry, while the minus 
sign applies to a ray moving towards the center. The 
constant h and the angle (f) between the direction of ob- 
servation and the direction defined by the observer and 
the center of symmetry are related by 



±\ I - 



(13) 



The angle cj) should be taken into account when one con- 
siders off-center observers (12] . 

In Ref.[13] it was shown that in the model ([5]) the 
inhomogeneity could mimic the dark energy and that 
^inh,o = 0.611q °o- It also emerged that the inhomo- 
geneity had dominated the universe quite recently, so it 
influenced only slightly the Doppler peaks and did not 
influence big-bang nucleosynthesis at all. Models of type 
I have been studied in Ref . [l5| , where it was shown they 
could be fitted to CMB data. 

Recently, lots of interest was attracted by the effect 
of redshift drift in cosmological models - the effect first 
noticed by Sandage and later explored by Loeb [l^ . The 
idea is to collect data from the two light cones separated 
by 10-20 years. It has recently been investigated for the 
LTB models [ij, backreaction timescape cosmology [l7| . 
and very recently for the axially symmetric Szekeres mod- 
els [l3| . Here we will consider this effect for the Stephani 
models. 

In order to do that we assume that the source does 
not posses any peculiar velocity, so that it maintains a 
fixed comoving coordinate dr = 0. The light emitted by 
the source at two different times te and te -I- St^ will be 
observed at to and to + 5to related by 



dt 



t„+st„ 



dt 



For small St^ and Sto we have 

Ste Sto 



a{te) a{to) 



(14) 



(15) 



A general formula for redshift which is valid for any cos- 
mological model reads as [Tsj l 



1 + z 



{Uak'')o 



(16) 



The redshift drift is defined as [l5j 

^ {Uak''){re,te+5te) _ (Mafc°)(re,te) 

{uak'^){ro,to + 5to) {uak<^){rQ,ta) 
which can be calculated using the expansions 

"9(Uafc") 



(17) 



(itaA:°)o = (uafc°)(r-o,io) + 

(Uafc")e = [Uak°-){re,te) + 

From (Uni) and ((H]) we have 
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(18) 



Using ([T5|) . and P^ . one can calculate the redshift 
drift p7)) for the Stephani universes as 
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(19) 



which with the help of ([8]) can be rewritten to the form 
Sz _ Ho 



1 + \H^o{^o - 1) 



X 







^r!oa-3(-+i) + (i-r!o)a-i- (1 + 2) 



(20) 



where a — a/ao and r = rao- 

Using all the above equations p5 |) - ([20|) we end up with 
the set of equations which combined together allow to 
find the redshift drift of any source at redshift z in the 
considered class of Stephani models II defined by the con- 
dition that k{t) — (3a{t) as 



Sz 
St 



4 + gg(r!o-l)5fg 

A + Hiino-iy 
(^^nod-'^^^+^'> + {i-no)d- 

4 + 1^4(1 + .), 



(21) 



(22) 



d7 



where r is the proper time (dr = dt/V) measured by an 
observer placed at fo and the last equation describes the 
propagation of the null geodesic. 

In the limit JIq ^ 1 =^ ^inh and w = (a flat 
Friedmann model filled with dust), the formula (|19|) re- 
duces to 



'St 



-Ho[il + zy^~{l + z)] , 



(24) 



which coincides with the standard Friedmann universe 
formula obtained by Sandage and Loeb [T^j. In 
Ref. 1211 the redshift drift as a function of redshift for 
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FIG. 1. The redshift drift flU^ as a function of redshift for the 
sperically symmetric inhomogeneous pressure Stephani model 
with ro = 0, w — 0, and Q,inh = 1 — f^o = 0.61. It shows 
that for large z the drift behaves alike in the void LTB mod- 
els. However, unlike in these void models, it becomes positive 
for small redshifts z G (0,0.34), reaching its highest value 
Ri 10"^^/year at « ~ 0.17. 



the ACDM model, the Dvali-Gabadadze-Porrati (DGP) 
brane model, the matter dominated model (CDM), and 
other three different LTB void models was presented. It 
has been shown that the redshift drift Sz/6t for ACDM 
and DGP models is positive up to z « 2 and becomes 
negative for larger redshifts, while it is always negative 
for LTB void models '22]. Using the formulas (CT. 
and ((23)) wc plot the redshift drift as a function of red- 
shift for the Stephani model with tq = Q, w ~ 0, and 
^inh = 1 — = 0.61 in Fig. [T] It emerges that for 
large redshifts the drift exhibits the behavior which is 
alike the redshift drift in the void models. However, 
unlike in the void models, the redshift drift becomes 
positive for z E (0, 0.34) and attains its highest value 
« 10^^^ /year for z ^ 0.17. This may allow to de- 



tect the difference between the spherically symmetric 
LTB models and spherically symmetric Stephani mod- 
els in future telescopes such as the European Extremely 
Large Telescope (E-ELT) (with its spectrograph CODEX 
(COsmic Dynamics Experiment)) [23], the Thirty Meter 
Telescope (TMT), the Giant Magellan Telescope (GMT), 
and especially, in gravitational wave interferometers DE- 
CIGO/BBO (DECi- hertz Interferometer Gravitational 
Wave Observatory/Big Bang Observer) Besides, it 
is worth mentioning that there is another test known as 
the cosmic parallax test [2l| which is due to anisotropic 
expansion and is strictly related to a non- vanishing shear. 
However, Stephani universes as shearfree should not ex- 
perience it - this may be another way to differentiate 
Stephani models with LTB models. 

Conclusion: In this paper we have discussed the 
Stephani models of pressure gradient shells which are 
complementary to the energy density varying shells of the 
Lemaitre-Tolman-Bondi models. The formula for red- 
shift drift Sz/St of any source at redshift z in the spe- 
cific class of Stephani models for both centrally and non- 
centrally placed observers has been obtained. We have 
shown that at least for the centrally placed observers, 
there is a subset of observationally viable Stephani mod- 
els which exhibit qualitatively different behavior of red- 
shift drift than the LTB void models. This difference 
may allow to test Stephani cosmology against LTB void 
cosmology in future experiments aimed to measure the 
redshift drift. 
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